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Abstract 

We predict a spatially localized magnetic domain wall oscillator upon the application of an 
external magnetic field and a DC electric current. The amplitude and frequency of the oscillator 
can be controlled by the field and/or the current. The resulting oscillator could be used as an 
effective microwave source for information storage application. 



1 



In a spin valve, a DC electric current generates a spin transfer torque which can control 
or alter the magnetization dynamics of the free layer. Above a critical current density, 
the spin transfer torque is able to switch the direction of the magnetization, generate spin 
wave excitations, and more interestingly, create a steady-state precessional motion of the 
magnetization of the free layer [l| . For a ferromagnetic metal, the spin transfer torque results 
in domain wall motion 0, 3], spin wave excitations j^, 5], and wall transformation from one 
type to another . What is yet to demonstrate is whether a DC current is able to create 
a spatially localized domain wall oscillator. A well controlled and spatially localized domain 
wall oscillator is very desirable for applications. For example, the oscillatory magnetic field 
from the stray fields of a localized domain wall oscillator can assist writing magnetic bits in 
recording media. Here, we show that a stable and localized domain wall oscillator is indeed 
possible by the combined applications of the magnetic field and the current. We determine 
the relevant parameters for this realization. 

Domain wall motion driven by a magnetic field has been well studied jsl, O, [h], 11]. 



When the magnetic field exceeds a critical value, the domain wall overcomes the pinning 
potential and begins to move to reduce the Zeeman energy. As a first approximation, the 
domain wall moves uniformly and the velocity of the wall is given by ^yHA /a js] where 
7 is the gyromagnetic ratio, H is the magnetic field along the wire, Ao is the wall width, 
and a is the damping parameter. When the magnetic field is further increased beyond 
the Walker-breakdown field, the wall motion is no more uniform; instead, the wall velocity 
becomes oscillatory [3]. However, the oscillation can not be sustained because the wall 
motion or wall oscillation continuously decreases the magnetic energy due to the presence of 
damping. After a typical time scale of nanoseconds, the domain wall either stops precessing 
or completely moves out of any finite regions, i.e., a spatially localized wall oscillation driven 
by a DC magnetic field alone is not realizable. 

We show below that the spatially localized wall oscillator is possible when a current is 
also applied. By properly choosing the direction and magnitude of the current density for 
a given magnetic field, the average velocity of the wall can be precisely controlled at zero, 
and thus a stable and spatially localized oscillator can be created. There is a key difference 
between the field-driven and current- driven domain wall dynamics: the energy damped in 
a period of wall oscillation can be compensated by the energy input from the spin transfer 
torque, but no compensation occurs for the magnetic field since the change of the Zeeman 
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energy is zero for a full cycle of oscillation. 

To determine the relevant parameters for the creation of the localized wall oscillator, we 
re dynamic equation of the magnetization in the presence of field and current 
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where H e // is the effective magnetic field including the external field, the anisotropy field, 
the magnetostatic field, and the exchange field; bj = Pje^B/ 'eM s and cj = where P is 
the spin polarization of the current; j e is the current density, hb is the Bohr magneton, e is 
the electron charge, and M s is the saturation magnetization. £ is a dimensionless constant 
which describes the degree of the nonadiabaticity between the spin of the non-equilibrium 
conduction electrons and local magnetization. 

We solve above dynamic equation in two ways. First, we analyze a simplified model based 
on the Walker's wall profile [8]: this enables us to analytically determine the condition for 
the formation of the wall oscillator. Numerical calculations are then followed to verify our 
analytical results. In Walker's model, the domain wall structure is characterized by two 
variables: the center position of the wall q(t) and the angle of the wall plane <fi{t). The wall 
width Ao is treated as a constant. With these simplifications, one can write Eq. (1) in terms 
of q(t) and <fi(t), 

A o + aq = 7 A # - cj (2) 
g-ctA o = 72vrM s A sin 20 - bj. (3) 

Eliminating q from above equations, we have 

(H -H w sm2<p) (4) 



r (1 + a 2 ) 

where we have defined field H = H+(a—^)bj/(jA Q ), and the Walker breakdown field H w = 
a2nM s . When H < H w , Eq.(4) has a steady state solution, i.e., = and sin 20 = H /H w ; 
this is the solution for the uniform motion of the wall with velocity: q = ( , yA H — Cj)/a. 
When Hq > H w , however, there is no steady-state solution; this is known as the Walker 
breakdown [8] . The direct integration of Eq. (4) yields 



tan <f>(t) = Cx + J 1 - C\ tan(cjt) (5) 



where d = H w /H and u = ly /H* - H 2 /{1 + a 2 ). Equation (5) indicates that the angle 
increases a it in one period = ixjuj. Thus the average angular velocity is 



$ = */t+ = (T^yV^o 2 -^- (6) 

By averaging Eq. (2) over one oscillation period and by using Eq. (6), we arrive at the 
condition for the localized oscillator (setting q — 0), 



The amplitude of the oscillation can be obtained by solving for q when q = 0. There are 
two solutions; their difference (divided by 2) is identified as the oscillation amplitude q amp - 
From Eqs.(2), (3) and (7), we have 

qamp ~ 2aH ' 1 j 

To further gain the insight on the solution of the localized wall oscillation, let us consider 
the change of wall energy in one period of oscillation. Multiply Eq (2) by q and multiply 
Eq. (3) by A o 0, and then substrate the resulting equations from each other, we have 

nA bj = a [ * dt(q 2 + A 2 <p 2 ). (9) 
Jo 

The right hand side represents the domain wall energy loss in one cycle, which has to be 
compensated by the work done by the spin torque on the left hand side. Note that the 
external field does not contribute any work in a complete cycle; this is the physical reason 
why the field alone is unable to sustain a localized domain wall oscillation. 

Up till now, we have considered the solution of localized wall oscillations in an ideally 
uniform film or wire. The center position of the wall oscillation is arbitrary as long as 
Eq. (7) is met. Any spatial variation of the parameters would lead to a drift of wall position. 
In order to stabilize the center of the wall oscillation at a desired location, one needs to 
design a structure that can suppress the drifting of the wall center but maintain the wall 
oscillation. At first, one might consider a local pinning to trap the oscillator, for example, 
by using a higher anisotropy material in a small region. However, we find that such local 
pinnings are not effective at all. If the pinning is strong, the wall oscillation is completely 
destroyed and a static domain wall will be formed at the pinning site. If the wall oscillation 



4 



persists over a weaker pinning potential, the wall center remains unstable against a small 
fluctuation of parameters. The reason is that the amplitude of the wall oscillation is several 
times larger than the wall width, see Eq. (8), so that the local pinning does not affect the 
oscillation significantly. We thus propose a scheme to stabilize the oscillation by designing a 
spatially varying damping parameter-this can be achieved via gradient doping of rare-earth 



impurities in ferromagnets 15| . We argue below that the wall oscillation is spatially stable 
in this design. 

Consider a spatial dependence of the damping parameter a(q). For a fixed current bj 
and field H, the center of the oscillation will be located at a certain position q = xo so that 
Eq. (7) is satisfied when a = az(xo). Then if there is a fluctuation, for example, the current 
density is slightly increased, the wall center will move along the direction of electron current. 
The equation of motion for the center of the wall is, from Eqs. (2) and (3), 

(2vrM s 



q ~ 7A # 



H 2 



da . . , 

m q (10) 



where H w <C H and a C 1 have been used. Clearly, if da/dq < 0, i.e., the damping 
parameter increases along the direction of the electron current (note dq < if b j is increased) , 
the drifting velocity of the center of oscillation q exponentially decays to zero. The wall 
oscillates around a new position near the original oscillation center, where Eq. (9) is still 
satisfied on average. 

Next we numerically solve Eq. (1) to validate our analytical results. We choose a magnetic 
wire whose width and thickness are sufficiently small so that the transverse wall is energet- 
ically favorable compared to the vortex wall. We also choose a = £ = 0.02 throughout the 
simulation; when a^(, there is no qualitative difference on the behavior of the wall motion 
except that the effective magnetic field has an additional term given by (a — £)6j/(7A ). In 
Fig. 1, we show the typical wall velocity q as a function of the magnetic field with and without 
the current. The linearly increasing of q at small fields represents the uniform steady-state 
motion of the wall. When Hq > H w , the average velocity decreases since the domain wall 
does a reciprocative motion, i.e., the wall motion is oscillatory but q is generally non-zero 
- the oscillator is not localized. At certain values of the field and current, we find q = 0, 
indicated by the mark "X" in Fig. 1. In the insert, we show the oscillation of the wall center 
position q(t) around a fixed point [x — 0). We notice that the current density required for 
the localized oscillator is relatively small, bj = 20m/ s or j e = 3.9 x 10 7 A/cm 2 (for P = 0.7), 



compared to the experiments on current-driven domain wall motion 0, 0, [l6,ll7|]; this makes 
experiments of searching for the localized oscillator easily accessible. 

We describe the stability of the localized oscillator by choosing a linearly varying damping 
parameter a(x) = 0.02 x (1 — x/L), where L is the wire length and —L/2 < x < L/2. We 
first consider the field H and the current bj such that the oscillator is localized at x = 0, see 
the point U X" in Fig. 1. When we slightly vary the current density or the magnetic field, 
the oscillator will relocate to a new center position x\ near x — 0. The time evolution of 
the displacement q(t) at this new position (not shown) is similar as the insets of Fig. (1). 
We show in Fig. 2 the position of the wall oscillation center, amplitude and frequency of the 
localized oscillator for the varying magnetic fields and currents. These results are in good 
agreements with the analytical results, Eqs. (6), (7), and (8). 

Finally, we emphasize that the localized domain wall oscillators proposed here is quite 
different from the previous work Q] where either an AC current or an AC magnetic 
field is used driving force. In those strong pinning potential via geometrical 

confinement is used to localize the domain wall and the oscillation of the wall is simply 
a response to the oscillatory external force (AC fields or currents). Our proposal here is 
to generate the wall oscillation by a DC magnetic field and to localize the oscillator by a 
DC electrical current. A spatial varying damping parameter can effectively stabilize the 
oscillator. 

This work is partially supported by DOE and Seagate Technologies Inc. 
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Figure Caption 



FIG.l (Color online) Average velocity of the domain wall v = q. The parameters are: 
M s = 800 (emu/cc), H K = 500 (Oe), A = 1.3 x 10~ 6 (erg/cm). Note we have chosen a 
reduced Breakdown field similar to the experimental value H w = 9 (Oe) The fitting 

curves are the analytical solutions of q from Eqs. (2) and (6), and the fitted wall width 
Ao = 20nm. Insets show the time evolution of the displacement q(t) at the point "X". 

FIG. 2 (Color online) Position of the wall oscillation center x\, amplitude q amp and frequency 
/ = u/tt of the oscillator as a function of the field and current: (a), (b) and (c) are for a fixed 
current, and (d), (e) and (f) for a fixed field. Damping parameter is a(x) = 0.02 x (1—x/L), 
where L = 0.8/xm. 
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